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Abstract 

It is shown that under lowering density or temperature a nucleon Fermi super- 
fluid can undergo a phase transition to a new superfluid state corresponding to 
<^ ' superposition of states with singlet-triplet (ST) and triplet-singlet (TS) pairing of 

^S) \ nucleons (in spin and isospin spaces) . Such states arise as a result of branching from 

O , one-gap solution of the self-consistent equations, describing ST pairing of nucleons. 

Q ■ The density and temperature dependence of the order parameters for new two-gap 

If} • solutions is determined in the model with Skyrme effective forces. 

<n : 
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As it is known, at sufficiently low temperatures a Fermi-liquid (FL) becomes unstable 
with respect to formation of Cooper pairs and passes into a superconducting (superfluid) 
state. Corresponding examples are an electron liquid in metals and superfluid phases of 
'"q ■ 3 He. There are some experimental facts in favour of superfluidity of atomic nucleus as 

well Since the study of superfluidity of atomic nuclei is difficult due to finiteness of their 

> . sizes, it is of interest to study superfluidity of infinite nuclear matter, especially in view of 
^ ■ astrophysical applications. This question has been considered in many works. However, 

earlier the normal-to-superfluid phase transition in nuclear matter has been studied only. 
In principle, at further lowering of temperature a superfluid FL of nucleons can again 
become unstable passing into a new superfluid state. Such states are characterized by not 
by one but a few order parameters and have a lower symmetry than the initial superfluid 
state. Thus, we are speaking of the superfluid-to-superffuid phase transitions in nuclear 
matter. We shall clarify the mechanism of appearance of new superfluid states as well as 
find the temperature dependence of the order parameters for these states. 

1. Basic equations. We study multi-gap superfluidity of nuclear matter using 
Landau's semiphenomenological concept of a Fermi-liquid. The basic formalism is laid 
out in more detail in Ref. @, where the phase transitions to one-gap superfluid states 
of symmetrical nuclear matter have been studied. The present work is an extension of 
Ref. and its main aspect is analysis of the phase transitions to multi-gap superfluid 
states. For simplicity we assume that the energy functional is invariant under rotations in 
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the configurational, spin and isospin spaces (for infinite uniform nuclear matter the spin- 
orbit interaction is equal to zero). Hence, superffuid phases (Cooper pairs) are classified 
by specifying the total spin of the pair, S = 0,1, the isospin T = 0, 1, their projections 
S z and T z on the z axis, and the orbital angular momentum L = 0, 1, 2, ... The possible 
values of the orbital angular momentum L for each value of S and T must be found 
according to the Pauli principle. Each superfluid phase is described by its own set of 
order parameters: at S = 0, T = by the scalar order parameter Aoo (singlet-singlet 
pairing of nucleons), at5 , = l,T = 0or5 = 0,T=lby the vector order parameter Ako 
or A fc, respectively {k = 1,2,3) (TS and ST pairing of nucleons) and at S = 1, T = 1 
by the tensor order parameter Ajfc (i,k = 1,2,3) (triplet-triplet pairing of nucleons). A 
nucleon superfluid FL is described by two distribution functions: the normal distribution 
function f KlK2 = Tr ga+a K1 and the anomalous distribution function g K1K3 = Tr ga K2 a K1 , 
where a+ and a K are the creation and annihilation operators of fermions with momentum p, 
spin (isotopic spin) projection er(r), k = (p, a, r), Trg... being the average of the operators 
with the density matrix of the system, g. The energy of the system is specified as a 
functional of the distribution functions / and g, E = E(f,g). It determines the fermion 
one-particle energy e and the matrix order parameter A of the system 

9E A n dE 

A ^ = (!) 



The matrix self-consistent equation for determining the distribution functions / and g 
follows from the minimum condition for the thermodynamic potential Q = — S + YqE + 
Y^aNa + Y^Nb (S is the entropy, N a , Nf, are the numbers of particles of species a and b, 
Y = 1/T, Y 4a = —fi a /T and Y^ = —fib/T are the Lagrange multipliers, \i a and //& are 
the chemical potentials of particles of species a and b, T is the temperature): 

/ = {exp(F ^ + %) + l} _1 = {exp(F + l}"' , (2) 




A+ -e 

Here the quantities e, A, I4 are, in turn, matrices in the space of the k variables, with 
^4ki« 2 = Y4 Tl 5 KlK2 (n = a, b) and the tilde stands for the transposition operation. Using 
the procedure of block diagonalization 0, |3| in an extended space 2x2, where the operator 
/ acts, one can express evidently the distribution functions / and g in terms of the 
quantities e and A. In what follows, we consider the case of symmetrical nuclear matter 
(equal numbers of protons and neutrons, or Y ia = Y ib = Y±, /i a — pb = A*)- in this case 
we obtain 

/ = lfl_| tanh MV 5 = _J_tanliM. A ; E = Jp + AA+,Z = e + l£ (3) 
j 2 1 E 2 l' y 2E 2 vs s Y Q K 1 
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We emphasize that the quantities f,g,£ and A in Eq.([|) are matrices in the spin, isospin 
and momentum spaces. Further we shall study the unitary states of superfluid nuclear 
matter, for which the product AA + is proportional to the product of the unit matrices 
in the spin and isospin spaces, AA + oc ct t . We consider as a particular case two-gap 
unitary states, for which the order parameters possess by the same symmetry properties 
on momentum: 

A(p) = A m (p)a 3 a 2 T 2 + A 03 (p)a 2 T 3 r 2 , A 30 (-p) = A 30 (p), A 03 (-p) = A 03 (p) (4) 

In this case the wave function of a Cooper pair describes the superposition of states with 
the triplet-singlet and singlet-triplet pairing of nucleons (TS-ST states): 

g(p) = 93o(p)^2T 2 + g 03 (p)a 2 r 3 T 2 , 03o(-p) = 9-io(p), gos(-p) = flte(p), (5) 

with the projections of total spin and isospin S z = T z = 0. 

2. Self-consistent equations. The case T=0. To derive the self-consistent equa- 
tions of a nucleon superfluid FL it is necessary to specify the energy functional, which we 
set in the form 

E(f,g) = E (f)+E int (g), E (f) = ± 5>(p)/oo(£), e (p) = |1 (6) 

V - 2m 
v 
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Eint(g) = yJ2i93o(P) V i(PiV)93o(q) +9h{p)V 2 {p,q)g 03 {q)} 
p,q 

Here m is the effective nucleon mass, foo is the coefficient of the product roUo in expansion 
of the distribution function / in the Pauli matrices, V\, V 2 are the anomalous FL interaction 
amplitudes, which have the symmetry properties Vi(—p, q) = Vi(p, —q) = Vi(p,q), i = 1, 2. 
The quantity m contains account of the normal FL effects and represents itself the mass 
of a free nucleon, renormalized by the normal FL interaction (in the Skyrme model the 
quantity m is given by Eq. (|10D). With allowance for Eqs. ([[]), (]|), we obtain expressions 
for the quantities A 30 , A 3 as functions of the quantities g 3 o and go 3 : 

AaoCp) = ^ E Vi& Q)9so(q), A 03 (p) = I £ V 2 (p, q)g 03 (q) (7) 

9 9 

Conversely, expressions for the distribution functions g^ 3 and g 3 Q as functions of the quan- 
tities £, A 30 , A 3 can be found from Eqs. (H), (|). With account of them, we get the 
self-consistent equations for determining the order parameters of two-gap unitary states 
of superfluid nuclear matter: 

Aao(p) = ~yz\ ^(p, q) |^ tanh + ^ tanh j , (8) 
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Here A± = A 30 ± A 03 , E± = y£ 2 + |A±| 2 . According to Eq. (J8|), considering a nucleon 
superfluid FL is characterized by two types of fermion excitations with the gaps A± in 
the spectrum. Eqs. (§) generalize equations of the BCS theory and contain the one-gap 
solutions with A 30 7^ 0, A 3 = (TS pairing) and A 30 = 0, A 03 7^ (ST pairing) as 
particular cases. 

Let us give an analysis of Eqs. (§), using the model representations of the BCS theory. 
Precisely, we assume that the interaction amplitudes V\ and V 2 are not equal to zero only 
in a narrow layer near the Fermi surface: |£| < 9, 9 ^ ep (in numerical calculations we 
set 9 = O.lep). Besides, we choose the effective Skyrme forces || as the amplitudes of NN 
interaction. Relation between the anomalous FL amplitudes V\, V 2 and the amplitudes of 
NN interaction has been set in Ref. 0. Taking into account it, for the quantities V\, V 2 
we have: 

ViM 0) = to(l ± x ) + \hn a {\ ± x 3 ) + ± Xl )(p 2 + g 2 ), (9) 

where n is the density of symmetrical nuclear matter, ti,Xi,a are some phenomenological 
parameters. Note that the amplitudes V±, V 2 contain no dependence from the parameters 
t 2 ,x 2 , because the amplitudes Vi,V 2 are the even functions of the arguments p, q (see 
details in Ref. M). There are sets of parameters U, Xi, a, which differ for various versions 
of the Skyrme forces (we shall use the Ska, SkM, SkM* and RATP potentials || as well 
as the SkP potential ||). In the Skyrme model the effective nucleon mass is given by the 
expression 

^ = £ + i( 3t ' + ^ 5 + 4 ^ < 10 > 

mo being the mass of a free nucleon. 

Let us consider the case T = 0. As a result of the above assumptions, we arrive at 
equations for determining the quantities A 30 = A 30 Qo = Pf), A 03 = A Q3 (p = p F ): 

A30 = 9i f_ 9 dt{(jr + -^)A 30 + ( JL - -^)A 03 }, (11) 
A 03 = 92 f e dt{{±- + -^-)A 03 + ( J- - A 30 } 



Here g>i j2 = : — and v F = 9 ^ 3 being the density of states at the 



vfVi, 2 (p = PF,q = Pf) , mp F 

and vp = — —^5 

4 2n 2 h 3 
Fermi surface of a proton (neutron) with the given spin projection. Eqs. (|ll|) allow 

to determine the order parameters A 3 o, Ao 3 in the TS-ST states of paired nucleons. In 

general case this can be done only numerically. Analytical consideration is possible if the 
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conditions A 30 , A 3 <C 9 are fulfilled (logarithmic approximation). Introducing the ratio 
of the order parameters a = A 30 /A 3 and performing integration in Eqs. ([11]) in this 
approximation, we arrive at the following equations for the quantities a and A03: 

1 — a 



1 , , 20 .1, 

— = 4 In — + 2(-ln 

9i A 03 a 

1 26 

— = 4 In- h 2 (a In 

#2 A 03 



1 + a 
1 — a 



1 + a 

Excluding A 03 from Eqs. fll2D , we obtain the equation 

1 + a 



-ln|l-a 2 |), (12) 
-ln|l-« 2 |). 



(p(oc) = , <p(oe) = 2(a ) In 

9i 92 a 



a: 



(13) 



Since (p m i n = <p(0) = —4 (the point a = is the point of removable discontinuity) and 
(Pmax = y?(±oo) = 4, then Eq. (|i~3|) has the solution for a, if the coupling constants <?i and 
(72 satisfy the unequalities —4 < l/pj — l/g 2 < 4. These unequalities have to be fulfilled 
in the logarithmic approximation for existence of the TS-ST mixed state. In this case, 
they impose certain restrictions on the coupling constants g± and g%. The sense of these 
restrictions is that the coupling constants in TS and ST pairing channels must be of the 
same order. Clearly that similar restrictions exist in a general case, when the conditions 
of the logarithmic approximation are not fulfilled. 

The results of numerical determination of the order parameters A 3 o(n), Ao 3 (n) as 
functions of density on the base of Eqs. (O) are displayed in Fig. 1. It is seen that at 
some critical density there appear two pair (±A 30 (n), A 03 (n)) of TS-ST solutions which 
thermodynamically are indistinguishable. We conclude that two-gap TS-ST superfluid 
states can arise in nuclear matter as a result of the density phase transition (at fixed 
temperature) from one-gap ST superfluid state (in the model with the Skyrme effective 
forces). Note that Eqs. ( |TT|) have two-gap solutions with A 30 ^ 0, A 03 ^ in the case 
of the SkM, SkM* and Ska potentials but have no such solutions for the RATP and SkP 
potentials. For the latter potentials the TS coupling constant is considerably larger than 
the ST coupling constant and, hence, the conditions for existence of two-gap solutions are 
broken. Concerning the potential SkP we remark that it has been used for the description 
of the pairing correlations in highly asymmetric nuclear matter || unlike to our case where 
we study superfluidity of symmetric nuclear matter. 

3. Temperature behaviour of the order parameters. Given in the section 2 
analysis relates to the case T = 0. It is clear, that if TS-ST states exist at T = then 
such states arise first at some critical temperature T tss t- To describe the temperature 
behaviour of the order parameters and the mechanism of appearance of the new solutions 
it is convenient to pass to the new independent variables x = A 30 + A 03 , y = A 30 — A 03 . 
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With allowance for Eq. (J8|), where we have replaced summation by integration, the self- 
consistent equations in terms of the variables x and y are written in the form 

^ = x\(x, T) + y\(y, T), = x\(x, T) - y\(y, T), (14) 

di E 



Kx, T) = J -| tanh ^,E= ^ + * 2 



Excluding the variable y from Eq. (|14]), for determining the function x = x(T) we obtain 
the equation 

d(x,T)-d(x-d(x,T),T) = D(x,T) = l, d{x,T) = 4 9i92X(x,T) - 9l - g 2 

92 - gi 

It is not difficult to see that the variables x and y are related by the formula y = xd(x, T). 
Therefore, the order parameters A 30 (T), A 3(T) can be found by the formulae: 

A 30 = \x{l + d(x, T)), A 03 = \x{\ - d(x, T)) (16) 

One-gap solutions are obtained from Eqs. (0) as solutions of the equations d(x, T) = 1 
(triplet-singlet) and d(x,T) = — 1 (singlet-triplet) while corresponding critical tempera- 
tures T ts , Tgt are determined from the equations d(0,T ts ) = 1, d(0,T st ) = —1. As follows 
from evaluations with the Skyrme forces, for all densities, where the superfTuid states exist, 
it holds g\ > Analyzing the behaviour of the function D(x,T), one can conclude that 
at temperatures T > T ts there is no one-gap solutions; at temperatures T st < T < T ts 
there is only one TS solution; at temperatures T tsst < T < T st the system is character- 
ized by two one-gap (TS and ST) solutions. Finally, at temperatures T < T tss t we have, 
besides the previous solutions, two new TS-ST solutions as well. 

To determine the critical temperature T tsst , at which TS-ST solutions arise first we 
use the following considerations. Clearly, TS-ST solutions arise as a result of branching 
from TS or ST branch of solutions Eq. (|15|) . In moment, when TS-ST solutions arise, the 
derivative D' x (x, T) vanishes in the critical point (xt ss t, T tS st)- If branching occurs from TS 
solution, i.e., d(x, T) = 1, then D' x = d x (2 + xd' x ) > for x > (since the function D(x, T) 
is even with respect to x, we can limit ourselves by the positive x; in this case d' x > 0) 
and the equation D' x = has no solution. If branching occurs from ST solution, i.e., 
d(x,T) = —1, then D' x = d' x (—2 + xd' x ) and the critical point (x tss t, Tt ss t) is determined 
from the equations 

d(x,T) = -1, xd' x (x,T) = 2 (17) 

Calculation of the second derivative D xx (x, T) in the critical point gives D xx (xt ss t, T tsst ) = 
0, i.e., the mechanism of branching of TS-ST solutions is formation of inflection on the 
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curve z = D(x,T tsst ) for x = x tss t- All said above allows to approximate the function 
D(x, t) for T < T st in vicinity of the inflection point in the form 



D{x, T) = A(x - x st (T)) 3 + C{T){x - x st {T)) + 1 (18) 

(x st (T) is the singlet-triplet branch of the equation D(x,T) = 1), where for the coeffi- 
cients A = ^D xxx (x tsst ,T tsst ), C{T) = D x (x,T)\ x=Xgt (r) it is not difficult to obtain the 
expressions 

A = - «) 3 - \d' x d'L ~ \<€ x )\ C(T) = -d' x (x st (T),T)(2-x st (T)d' x (x st (T),T)). 

Retaining in expansion of the function C(T) the linear with respect to T tsst — T term, we 
obtain the formulae, describing the behaviour of TS-ST solutions near the critical point 



42(r) = x tsst ± yJp(T t8St -T)/A, 13 = 2d" xT - d' T (d' x + wC)l (19) 

In this case, it is necessary that (3/ A > 0. 

The results of numerical determination of the order parameters A 30 (T), A 03 (T) on 
the base of Eqs. (|T5|), ( p!6|) are displayed in Fig. 2. It is seen that two pairs of TS- 
ST solutions (±A 30 (T), A 03 (T)) exist for the temperatures T < T tsst and besides, in the 
critical point, in accordance with Eq. (|16|), A 30 = 0, A 03 = Ag 3 (Aq 3 being the one-gap 
ST solution). We recall, that we consider the case x > 0. Existence of two pairs of the 



solutions (±A 30 (T), A 03 (T)) follows from the fact that, if Eqs. (|14]) have as a solution the 
pair of functions (x(T),y(T)), then after substitution x — > — y, y — > — x we obtain the 
new pair that is a solution too. In the case x < we would be able to obtain another two 
pairs of TS-ST branches (±A 30 (T), — A 03 (T)). 

Conclusion. Thus, we pointed out the possibility of arising of two-gap unitary states 
in superfluid nuclear matter corresponding to the superposition of ST and TS pairing of 
nucleons. The self-consistent equations for these states differ essentially from the equations 
of the BCS theory and contain the one-gap solutions (ST and TS) as some particular cases. 
Analysis of the self-consistent equations at the temperature T = in the logarithmic 
approximation shows that new states can arise only under quite specific restrictions on the 
coupling constants, which describe interaction of nucleons in ST and TS pairing channels. 
Since the constants of the effective interaction depend from density, there exist such 
density intervals for which either one-gap or two-gap superfluidity of nuclear matter is 
realized. Calculations with the Skyrme interaction being chosen as a potential of NN 
interaction indicate that two-gap unitary states arise in nuclear matter as a result of the 
phase transition in density from one-gap ST superfluid state (results have been obtained 
for the Ska, SkM and SkM* potentials). Analysis of the temperature behaviour of the 
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order parameters shows that new two-gap states can arise under lowering temperature 
from the superfluid state with ST pairing of nucleons as well. Thus, we conclude that a 
superfluid nucleon FL under lowering density or temperature can undergo the new phase 
transition leading to emergence of two-gap superfluid states. Among the other problems 
we note here the study of influence of asymmetry |J and real bound states (deuterons) 
[0. U on the thermodynamic properties of multi-gap superfluid states in nuclear matter. 
The authors thank BMBF for partial financial support of this work. 
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Figure captions: 

Fig. 1. Order parameters A 30 , A 3 vs density of symmetric nuclear matter at T = 
in the case of TS (A 30 ^ 0, A 03 = 0), ST (A 30 = 0, A 03 ^ 0) and mixed TS-ST (A 30 ^ 
0, A 03 7^ 0) states of paired nucleons for the SkM* (a), SkM (b), Ska (c) Skyrme forces; 
tsst(ts + ), tsst(ts-) and tsst(st±) are notations for the dependences of the TS and ST 
order parameters A 30 and A 03 in two pair (±A 30 (n), A 03 (n)) of TS-ST solutions of Eq. 
(|15"D. Branching of TS-ST solutions occurs at n = 0.049 fm -3 for the SkM* potential, at 
n = 0.052 fm -3 for the SkM potential and at n = 0.064 fm -3 for the Ska potential. 

Fig. 2. Order parameters A 3 o, Ao 3 vs temperature. Notations are the same as 
in Fig. 1. Calculations have been performed for the SkM* (a), SkM (b) potentials at 
n = 0.04 fm -3 , for the Ska (c) potential at n = 0.05 fm -3 . Branching of TS-ST solutions 
occurs at T = 0.162 MeV for the SkM* potential, at T = 0.179 MeV for the SkM potential 
and at T = 0.203 MeV for the Ska potential. 
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